Cosmic-ray electron signatures of dark matter 
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There is evidence for an excess in cosmic-ray electrons at about 500 GeV energy, that may be 
related to dark-matter annihilation. I have calculated the expected electron contributions from a 
pulsar and from Kaluza-Klein dark matter, based on a realistic treatment of the electron propagation 
in the Galaxy. Both pulsars and dark-matter clumps are quasi-pointlike and few, and therefore their 
electron contributions at Earth generally have spectra that deviate from the average spectrum one 
would calculate for a smooth source distribution. I find that pulsars younger than about 10 5 years 
naturally cause a narrow peak at a few hundred GeV in the locally observed electron spectrum, 
similar to that observed. On the other hand, for a density n c = 10 kpc~ 3 of dark-matter clumps 
the sharp cut-off in the contribution from Kaluza-Klein particles is sometimes more pronounced, 
but often smoothed out and indistinguishable from a pulsar source, and therefore the spectral shape 
of the electron excess is insufficient to discriminate a dark-matter origin from more conventional 
astrophysical explanations. The amplitude of variations in the spectral feature caused by dark 
matter predominantly depends on the density of dark-matter clumps, which is not well known. 

PACS numbers: 95.35.+d; 96.50.sb 



I. INTRODUCTION 



Only about 1 per cent of galactic cosmic rays are elec- 
trons, but their properties are of particular interest be- 
cause they are very radiative at high energies and thus 
quickly loose their energy. The cosmic-ray electron spec- 
trum is therefore softer than that of cosmic-ray nuclc- 
ons, and high-energy electrons are few. For a long time 
only emulsion-chamber data of the electron flux above 
100 GeV were available which were well represented 
by a power law N(E) oc E~ 3 - 2 , but the energy resolu- 
tion and statistical accuracy were limited. Recently data 
obtained with the ATIC balloon experiment were pub- 
lished, which show an excess of galactic cosmic-ray elec- 
trons at energies between 300 and 700 GeV 0. At TeV 
energies the electron flux appears to drop off rapidly [l| . 
This excess is indicative of a previously unknown indi- 
vidual source of high-energy electrons, which could be a 
nearby supernova remnant [J] , a pulsar [f| , a microquasar 
Q , or an annihilation site of dark-matter particles of the 
Kaluza-Klein type [?J • 

The dark-matter interpretation is particularly appeal- 
ing, because the PAMELA collaboration has reported 
an increase in the cosmic-ray positron fraction above 
20 GeV, suggesting the existence of a local source of both 
positrons and electrons Q . Dark-matter annihilation re- 
sulting in electron-positron pairs is possible in a number 
of models d, Ell El, E2], and will produce a spectrum 
dominated by a delta functional at the mass of the dark- 
matter particle EII . Indeed, the ATIC team finds that a 
Kaluza-Klein particle with mass 620 GeV fits their elec- 
tron data just fine, when the expected electron source 
spectrum is propagated in the Galaxy using the GAL- 
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PROP code [l4[. Hall and Hooper suggest that high- 
precision measurements of the electron spectrum be con- 
ducted with atmospheric Cherenkov observatories such 
as VERITAS and HESS, that would permit a discrim- 
ination between the dark-matter and pulsar hypotheses 

m. 

If dark-matter annihilation is responsible for the ex- 
cess in 600-GeV electrons, then a substantial boost fac- 
tor is required to match the observed electron flux which 
requires that the dark matter be concentrated in dense 
clumps. The electrons would then be injected into the 
Galaxy only at the location of those clumps, which intro- 
duces substantial variations in the electron flux through- 
out the Galaxy and can significantly modify the ob- 
served electron spectrum, as was shown in similar stud- 
ies of electron propagation from supernova remnants 
EI E3, E3- The GALPROP code implicitely assumes 
a smooth source distribution on account of its using a 
finite-difference algorithm on a coarse grid, and therefore 
it will not properly describe those fluctuations. 

Here we study the propagation of relativistic electrons 
from localized sources. 



II. THE PROPAGATION OF RELATIVISTIC 
ELECTRONS 

The effects of the spatial structure of the electron 
sources appear only at higher particle energies, at which 
the radiative loss time is short. Therefore we may treat 
the propagation of electrons at energies above 50 GeV 
with a simplified transport equation, 

^-JL(bE 2 N)-DE a V 2 N = Q (1) 

with which we consider continous energy losses by syn- 
chrotron radiation and inverse Compton scattering, a dif- 
fusion coefficient D E a dependent on energy, and a source 
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term Q. Throughout this paper the propagation param- 
eters have the values 



1 

bE 



2.6 • 10 15 s, 



(2) 



corresponding to synchrotron losses in a 11 pG galactic 
magnetic field and Compton scattering of the CMB. The 
diffusion coefficient is chosen as required in cosmic-ray 
propagation without continuous reacceleration [l^, [2(| 
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Green's function for this problem is [2l| 

G = ^ exp 

bE 2 (4ttA) 3/2 

where is a stepfunction and 

D {E 11 - 1 - E' a ~ x 
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A = 
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(5) 



In the case of discrete sources the injection term Q is a 
sum over all such sources. For an individual source we 
can write 
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FIG. 1: Examples of electron spectra measured near Earth 
for different ages of the pulsar (cf. Eg. I10p . The thin dotted 
line denotes the pulsar source spectrum as in Eq.(8] The three 
thick lines give the particle spectra for three different values 
ages parameter f , but the same distance parameter, p. 
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and obtain the current (t = 0) contribution of that source 
to the electron density at distance r as 
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A. Electrons from pulsars 

The spectrum of electron escaping from a pulsar is not 
well known. A simple parametrization may be in order 
[IH that describes the differential production rate of elec- 
trons by a pulsar 



f p (E) = E-^ exp 



E 

E~r. 



g p (t)=@(t + r) 



(8) 



(9) 



where we use a stepfunction to account for the finite 
age, r, of the pulsar. Obviously, the spectrum of the ex- 
cess electrons depends chiefly on the scaling parameters 
i = bE c r and p = (r 2 [1 - a] 6)/(4£) J B^ 1 ), which com- 
pare the age, r, with the energy-loss timescale,l/(6£' c ), 
and the distance, r, with the diffusion length within one 



energy-loss time. At time t = and distance r we observe 
the differential density of electrons as 
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where C absorbs all constants and 

if #£>1, 



oc 



(11) 



In Fig. [T] we show possible electron spectra near Earth, 
that may result from a pulsar that produces electrons 
with spectrum ([5]), indicated by the thin dotted line. To 
be noted from the figure is the sharp peak in the elec- 
tron spectrum that is entirely a propagation effect. The 
distance to the pulsar is assumed as r — 700 pc; for our 
propagation parameters the distance parameter is then 
p = 0.23, meaning electrons observed at 600 GeV can 
well reach Earth within one energy-loss time, whereas 
electrons beyond a few TeV can not. 

If £ = 1, the pulsar age is the same as the energy- 
loss time at 600 GeV, 140,000 years for our propaga- 
tion parameters. Electrons at lower energy may not have 
enough time to not reach Earth because their propaga- 
tion range scale oc E a , thus causing a sharp low-energy 
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cut-off in the observed spectrum. Given the local super- 
nova rate we expect about one pulsar born within 1 kpc 
and 10 5 years, so given the energy-loss time at 600 GeV 
one pulsar at r — 700 pc is realistic. However, the pul- 
sars need a few million years to propagate 1 kpc above 
or below the galactic plane, and so electron injection by 
the pulsar must taper off after about 10 5 years to avoid 
a high flux of 100-GeV electrons. 



III. ELECTRONS FROM DARK MATTER 

Having established in the preceding section that for re- 
alistic parameters a pulsar can produce a narrow peak in 
the local electron flux, I will now discuss electron spec- 
tra that may arise from dark-matter annihilation. As 
substantial boosting factors of a few hundred are needed 
[22| . the dark matter is most likely organized in a num- 
ber of individual high-density clumps. The size of the 
clumps is not relevant for us, as long as it is much smaller 
than the propagation range of 600-GeV electrons, a few 
hundred parsec. Likewise, the large-scale distribution of 
the clumps in the Galaxy does not matter, unless the 
clump density varies on scales similar to the electron 
propagation range. I will therefore assume the clumps 
to be randomly distributed in space with constant den- 
sity n c . Dark-matter annihilation should proceed at a 
constant rate and can, in the case of Kaluza-Klein par- 
ticles, produce electrons with a source spectrum that is 
dominated by a delta-functional at the particle mass. For 
each clump, the differential source rate of electrons can 
then be described by the functions (cf. equation©) 
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(12) 



The electron spectrum observed at distance r from the 
clump is then 



N dm = 



CO(E c -E) 
E 2 A 3 / 2 



eXP, ~4A 



(13) 



where C absorbs the constants, O is a stepfunction, and 
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(14) 



The total electron spectrum is then obtained by summing 
the contributions from all clumps. If the clump density 
n c is high, the dark-matter distribution is effectively ho- 
mogenous. Then the total electron spectrum is governed 
by a cooling tail. 



N dn 



'dr Airr 2 N dm = C' Q{E ' 
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E 2 



(15) 



This is the case implicitely (and tacitly) assumed when 
using the standard GALPROP code, and it is presented 
in many publications [e.g.[2|. 
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FIG. 2: The range of electron spectra measured near Earth 
for different clump densities of dark matter, added onto the 
galactic electron background with spectrum oc E~ A ' 2 . The 
dotted line indicates the mean spectrum (cf. Eq. 1150 , The 
dash-dotted line denotes a randomly selected spectrum as an 
example of what may be observed. The shaded areas indicates 
the range in which we find the electron flux in 68% (dark 
gray), 90% (medium gray), and 99% (light gray) of all cases. 



The question arises at what density the spectra start 
deviating from that for the homogeneous case (Eq fT5l) 
and what the observed spectral shape might be. 

I have randomly placed in the Galaxy dark-matter 
clumps with constant density n c and summed their elec- 
tron contribution according to equation (|13|) . using E c = 
600 GeV. The resulting electron spectrum from dark mat- 
ter is added to the generic galactic electron flux, for which 
a spectrum oc E~ 32 is assumed. On average, the dark- 
matter component at 600 GeV has twice the flux of the 
galactic electron background. 

Figure [2] shows the resulting total electron spectra for 
three different clump densities. In all three panels the 
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dotted line indicates the mean spectrum according to 
Eq. (fl~5|) , and the dash-dotted line denotes a randomly 
selected spectrum out of the 5000 that were calculated. 
To be noted from the figure is that the dark-matter hump 
often doesn't look different from the hump a pulsar would 
produce, ft can be fairly roundish and lack the sharp cut- 
off at E c . The reason is that electrons at an energy very 
close to E c must be very young, because they haven't lost 
a significant fraction of their energy, and can therefore 
only come from a very close dark-matter clump. Even 
for n c = fOO kpc~ 3 we expect only one clump within 
1 40 parsec, and thus electrons at E ~ E c may not reach 
us. On the other hand, a very close dark-matter clump 
would produce a dominant spike at E ~ E c . 

The shaded areas indicate the range of flux for three 
different probabilities. In 68% of all cases the electron 
flux is within the dark gray region, the medium gray area 
corresponds to 90% and the light gray to 99% probability. 

IV. SUMMARY AND DISCUSSION 

The ATIC collaboration has measured an excess in 
cosmic-ray electrons at about 500 GeV energy Q , which 
may be related to dark-matter annihilation. In this pa- 
per I have calculated the expected electron contributions 
from a pulsar and Kaluza-Klein dark matter. My empha- 
sis is on a realistic treatment of the electron propagation 
in the Galaxy, for which I use analytical solutions to the 
electron transport equation. The commonly employed 
GALPROP code implicitely assumes a smooth distri- 
bution of the electron sources, because it uses a finite- 
difference algorithm on a grid. 

The findings can be summarized as follows: 

• Pulsars younger than about 10 5 years naturally 
cause a narrow peak at a few hundred GeV in the 
locally observed electron spectrum. A single pul- 
sar could therefore explain both the electron ex- 
cess measured with ATIC and a similar excess in 
positrons, evidence for which at 50 to 100 GeV was 
obtained by the PAMELA experiment 0. The 
pulsar hypothesis does require that pulsars with 
ages 10 5 to 10 6 years leak significantly fewer elec- 
tron/positron pairs, otherwise they would provide 
a very strong contribution in the 50 to 300 GeV 
band that is not observed. 

• Dark-matter annihilation occuring predominantly 



in dense clumps will produce a feature in the 
local electron spectrum, that deviates from that 
expected if the dark matter were smoothly dis- 
tributed. The sharp cut-off in the contribution 
from Kaluza-Klein dark matter is often smoothed 
out, and the spectral feature would be indistin- 
guishable from a pulsar source, even if the energy 
resolution of the electron detector were perfect. 
The spectral shape of the electron excess is insuf- 
ficient to discriminate a dark-matter origin from 
more conventional astrophysical explanations, con- 
trary to a recent claim [151 ]. 

• While the mass of the dark-matter particle may be 
misestimated by only 20% or so, the amplitude of 
the electron excess can vary by more than a factor 
of 2 for a clump density n c — 10 kpc -3 , and the 
required boost factors will be misestimated by the 
same factor. 

• All variations in the amplitude and spectral shape 
of the dark-matter contribution to the local elec- 
tron flux depend on the density of dark-matter 
clumps, the variations being larger for smaller 
clump densities. 

If the clump density is the decisive parameter determin- 
ing the amplitude of spectral variations in the electron 
excess from dark-matter annihilation, then it is of prime 
interest to estimate that number. Dark-matter clumpy- 
ness provides the boost factors required by the ATIC 
data, which may be further increased by Sommerfeld cor- 
rections [23l . [24j | . Simulations of structure formation in 
cold-dark-matter cosmologies show clumping on a variety 
of scales, but the boost factors are generally very mod- 
erate [25|, [2^| . The clump density in those simulations is 
generally smaller than assumed in this paper, but that 
may be due to limited numerical resolution. The density 
of simulated particles in Via Lactea II barely exceeds 
10 3 kpc -3 at the solar circle, and therefore the clump 
density will unavoidably be much lower. In any case, 
if the clump density is indeed significantly lower than 
n c = 10 kpc -3 , then the dark-matter scenario can also 
not be distinguished from a pulsar origin by studying the 
high-latitude diffuse gamma-ray emission from the excess 
electrons, because the gamma-ray intensity distribution 
in the dark-matter scenario becomes similarly patchy as 
in the pulsar case. 
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